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Abstract 

We estimate the lattice artifacts in loop correction perturbatively for domain- 
wall QCD with infinite number of extra flavors. We find that there appear no 
0{a) errors in renormalization factors of quark wave function, mass and quark 



^^ \ bilinear operators at one and two loop level with off-shell quark momentum. 

f"^ \ Our proof is based on even or oddness of the quantum correction in terms of 

cn ' the quark external momentum and mass, and it can be extended to any loop 

y^ ■ level. 
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I. INTRODUCTION 

Syinanzik's improvement program |jll| applied to on-shell quantities p| attempts to elim- 
inate cut-off dependence order by order by an expansion in powers of the lattice spacing a. 
In order to improve the Wilson quark action to 0{a) in lattice QCD, this requires adding the 
0{a) "clover" term P]. Quark operators also have to be modified by 0{a) counter terms, 
which generally involve new operators of higher dimension 0-|^. In perturbation theory, the 
tree-level value of the clover coefficient and those of the counter terms of quark operators are 
sufficient to remove terms of 0{g^aloga) in on-shell Green's functions evaluated at one-loop 
level [0. However there still remains 0{g'^a) terms and new counter terms are needed to 
ehminate them |]3,||]. 

The domain-wall fermion formulation was originally proposed by Kaplan to define 
lattice chiral gauge theories with the introduction of many heavy regulator fields as an 
extension of the Wilson fermion. This original idea and the later work about the Chern- 
Simons currents |]TD[ is further developed by Shamir [n|,|T2| into a simpler form and is applied 



to lattice QCD (DWQCD) anticipating superior features over other quark formulations: no 
need of the fine tuning to realize the chiral limit and no restriction for the number of flavors. 
Recent simulation results seem to support the former feature non-perturbatively [^^-[T6|. It 



is also perturbatively shown that the massless mode at the tree level still remains stable 



against the quantum correction [|22|. In spite of the existence of the Wilson term in the 
action it does not directly affect the mass term; there is no additive mass correction. 

Related to this good chiral property and disappearance of additive mass correction from 
the Wilson term the DWQCD is expected to have no next leading 0{a) errors. It was argued 
in Ref. |T^ with an intuitive discussion of good chirality that the leading discretization errors 
are 0{a'^) and there is no 0{a) errors. They also argue that the scaling behavior of numerical 
results for Bk H^ and strange quark mass IQ are consistent with this nature. For the two 



dimensional Gross-Neveu model it is reported in Ref. [|I^ that the physical observables do 

not have 0{a) errors by a nonperturbative analysis of the effective potential at large A^ limit. 

Although the lattice axial symmetry was discovered |T^ for the Dirac operator which 

satisfy the Ginsparg- Wilson relation |[T9|j2( 



{D,^^} = aDR^^D (I.l) 



and it was shown to exist in the effective theory of DWQCD |21||, disappearance of 0{a) 
term is not necessarily trivial. Because the lattice axial transformation of Liischer type 
contains 0{a) term within the transformation itself, 

5i, = ^^{l-aRD)iP, # = ^75 (1-2) 

and the operators with exact chirality is given by using the projection 

Ot = iiT{l- '^RDJ ij. (1.3) 

This projected operator can be written in terms of the boundary fermion field of the DWQCD 
Or = q^q [0 , where the boundary fields are given by 

q=(l- ^Rd] ^, q = ^ (1.4) 



and are transformed under the axial transformation ( p2D as 

6q = 75g, 6q = q-f5- (1-5) 

In the context of the DWQCD, where the exact axial symmetry does not exist because of 
infinite number of "unphysical" fermion fields, the chirality of the operator is understood 
by the axial Ward-Takahashi identity given by Furman and Shamir with an appropriately 
defined axial transformation [|l^, 



V, {A;{x)0) = 2m {Jt{x)0) + {Jt,{x)0) + ^ {S^O) , (1.6) 

where A" is an axial current, Jg is a pseudo scalar density and J^ is an explicit breaking 



term. It was shown in Ref. [12] that the exact axial Ward-Takahashi identity is satisfied for 



boundary quark operators in the infinite flavor limit, 

V^ (A^(x)O(g, g)) = 2m {J^{x)0{q,q)) + i {dlOiq^q)) . (1.7) 

The above operator Or = q^q transforms in the same way as in the continuum under 
the transformation and the exact chirality can be given in order to keep the identity. Al- 
though the axial Ward-Takahashi identity is shown to be realized nonperturbatively for a 
well behaved gauge field configuration |12[, it is still interesting to see the perturbative un- 



derstanding of this good chirality, especially the disappearance of the 0{a) terms since we 
can see the operator mixing structure definitely. 

In this paper, we estimate the lattice artifacts in loop correction perturbatively for 
DWQCD with infinite width of fifth dimensional direction. We show that there appear 
no 0{a) errors in renormalization factors of quark wave function, mass and quark bilinear 
operators at one and two loop level. Although our instrument is perturbation theory our 
proof is based on even or oddness of the quantum correction in terms of dimensionful quan- 
tity such as the quark external momentum and mass, and it can be generally extended to any 
loop level. We notice that we need not to set the external quark momentum to the on-shell 
value to eliminate the 0{a) errors as in the case of the Wilson fermion with clover term. The 
0{a) errors automatically vanish for off-shell quarks like naively discretized lattice fermion. 

This paper is organized as follows. In Sec. |I| we introduce the DWQCD action and 
the Feynman rules relevant for the present calculation. In Sec. |T| we estimate the lattice 
artifacts for the quark self energy by expanding the correction in terms of the external 
quark momentum and mass. We start our calculation at one loop level. We notice that the 
quantum correction can be classified into odd function of quark mass and momentum. Since 
the relevant term for renormalization is given as a leading order and the 0{a) error is a next 
to leading order term in expansion, we can see the absence of 0{a) term quite easily. This 
nature of even or oddness is not specific to one loop level we can apply our procedure to 
the two loop and any loop level diagrams. Sec. ^ is devoted to the estimation of the quark 
bilinear operator effective vertex at one and two loop level. Our conclusion is summarized 
in Sec. |V|. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(A^c) gauge group with the gauge coupling g. We set 
number N of the regulator field or the length of fifth dimensional direction in domain-wall 
fermion to infinity. 



II. ACTION AND FEYNMAN RULES 



We adopt the Shamir's domain- wall fermion action |TT|, 



N 



5l 



DW 



EE 

n s=l 



- J2 {'^{n)s{-r + -rf,)Uf,{n)-^{n + /i), + ^(ra),(-r - -f^,)Ul{n - ^)i){n - /x), 



n)s{l + 75)V'Hs+i + i^{n)s{l - 75)^(^)s-i) + (M - 1 + 4r)^(n),V'(^) 



+ m^ (V'(n)7vPR?/'(n)i + ^{n)iPL^{n)N 



(ii.i: 



where n is a four dimensional space-time coordinate and s is an extra fifth dimensional or 
"flavor" index, the Dirac "mass" M is a parameter of the theory which we set < M < 2 
to realize the massless fermion at tree level, m is a physical quark mass, and the Wilson 
parameter is set to r = —1. It is important to notice that we have boundaries for the flavor 
space; 1 < s < A^. In our one-loop calculation we will take N -^ oo limit and our proof of 
vanishing 0{a) error is valid only in this limit. Pr/l is a projection matrix 



R/L 



1±75 



(11.2) 



For the gauge part we employ a standard four dimensional Wilson plaquette action and 
assume no gauge interaction along the fifth dimension. 

In the DWQCD the zero mode of domain-wall fermion is extracted by the "physical" 
quark field defined by the boundary fermions 



q{n) = PRi){n)i + PLi'{n)N, 
q{n) = ip{n)NPR + i){n)iPL. 



(11.3) 



We will consider the QCD operators constructed from this quark fields and we calculate the 
lattice artifacts in the Green functions consisting of only the "physical" quark fields. 

Weak coupling perturbation theory is developed by writing the link variable as U^,^ = 
g«9^M{^+A/2) g^j^^ expanding it in terms of gauge coupling g. The free domain- wall Dirac 
operator is given as a leading term 



Pj>)st = X1^7Msinp;,5st + (W'^ip) 



where the mass matrix is 



w+{p) = {w-{p)y 

M+ = {M-f = 
W{p) = 1-M 



rriM^ 



-W{p) 



s,t 



PR+(W-ip)+mM-) Pl 



(11.4) 



-W{p) 



-W{P)_ 



rE(i 



cos 



P/.)- 



(11.5) 

(11.6) 
(11.7) 



The domain-wall fermion propagator is given by inverting the Dirac operator (p.4| ) for 

A^ ^ oo limitQ 



Spip, m)st = -ilf, sinpf, + W + niM Gr{p, m)utPR 
-27. sin p. + W^ + niM^ Gl{p, m)utPL, 

^ su 

where sum over the same index is taken implicitly. Gr/l is given by 



(11.8) 



GR{p,m)st = GLip,m)N+i-s,N+i- 



A{1 



m 



m , 



F 

~a{N~s+t) 



1 - iye-")e-"('^-^-*) + (1 



We' 



-a{s+t) 



^ ^-a{N+s-t)^ ^ ^^-a\s-t\ 



cosh a 



l + W^ + E^sm^P^. 



2\W\ 

A= (2Vrsinha)-^ 
F = l-e'^W 



m^{\ 



iye-"i 



(11.9) 

(11.10) 

(11.11) 
(11.12) 



Note that the argument p of factors a and W is suppressed in the above formula. Since we are 
interested in the Green functions constructed with "physical" quark fields the above fermion 
propagator appears only as internal quark line. In order to construct the whole Green 
function we need other three types of fermion propagators which connects two "physical" 
quark fields and "physical" quark with fermion field of a general flavor index, 



-Z7^sinp^ + (1 - W^e" 



vn 



where 



^^(^)^(-^)^=-(l-e"W^) + m^(l-iye-) 
\q{p)i){-p)s) = rio{p,m)L{p)s + rie{p,m)R{p), 
(V^(p)sg(-p)) = R{p)s'no{p,m) + L{p)sr]eip,m) 



R(p)^ = e-"(p)(«-i)p^ + e-"(p)(^-^)Pi, 
r]o{p, m) = — (i-ff, sinp^ - m (l - 1^6^°)) 



= Sg{p,m), 



Ve{p,m) 



m 



(^7/. 



smp^ — m 



(1 



We- 



(11.13) 

(11.14) 
(11.15) 



(11.16) 
(11.17) 

(11.18) 
(11.19) 



Here we notice that Sq{p,m), rio{p,m) are odd function of {p,m) and L{p)s, R{p)s, Ve{p,^ 
are even function since W{p), a{p) are even function of momentum p^, 



^ Although the N ^ 00 Hmit should be taken after the loop momentum integral in principle, it 
is easy to see that the limit can be taken before the integral by neglecting terms proportional to 



-aN 



We take the form of free fermion propagator in the N ^ 00 limit from the beginning. 



Sq{-p, -m) = -Sg{p, m), rio{-p, -m) = -r]o{p, m), 
rie{-p, -m) = ?7e(p, m), L{-p)s = L{p),, R{-p)s = Rip)s- 



(11.20) 
(11.21) 



It is important for our later calculation that the quark-fermion propagator ( |II.14| ) and ([II.15|) 
can be classified into even and odd definite part with different "damping factor" L{p)s and 

Rip)s- 

In the perturbative calculation of loop diagrams we take external momenta and quark 
masses much smaller than the lattice cut-off, so that we can expand the external quark 
propagator in terms of quark momenta and masses. In this paper we adopt the following 
form of expansion in order to extract the one particle irreducible (IPE) vertex function from 



the loop diagram 23 , 




Sq{p,m) = 


^ ^° n 1 n(n^)) 




{q{p)i^{-p)s) = 


1 _ ,,,2 


tp + [1- WQ)m 


mp)sq{-p)) - 


= .Rip)s^eip,m)- L{p)s^oip,m)[ 



^o{p, m)R{p), 



Wa 



+ (1 



WQjm 



(11.22) 
(11.23) 
(11.24) 



where, Wo = W{0) = 1— M. We notice that the quark-fermion propagator (|II.23| ) and (|II.24|) 
are given as multiplication of the free quark propagator ( |II.22| ) with the linear combination 
of damping factors, whose coefficients ^^{Pi^i^oiP-,^ are given in an expanded form. 



ie{p,m) 
io{p,m) 



1 + 0{a^\ 
Wo . J 



Wn 



+ 0{a^). 



(11.25) 



We call this linear combination as external quark line factor in the following. An explicit 
form of the next to leading term is not important for our calculation. What is relevant for 
the later discussion is the fact that ^e(p, '"^) and ^o{p,fn) is an even and odd function in 
{p, m) and the even-odd structure of this external line factor is strictly characterized as a 
coefficient of damping factors L{p)s and R{p)s, say the even function ^g is accompanied with 
Lg and the odd function ^o with Rg in (p!I.23|) and opposite pairing in ( [II.24|) . We did not 
expand L{p) and R{p) because their next to leading terms do not affect the characteristic 
form of (|II.23|) (|II.24|) but only shift the value of ^e/o and the damping rate a. 

Since our gauge part is same as that of the usual Wilson plaquette action, the gluon 
propagator can be written as 



Gtip) 



1 



4 sin^ p/2 



Sfiu - (1 - a 



4sinp^/2sinpy/2 
4sin^p/2 



(11.26) 



where sin^p/2 = X)/^ sin^p^/2. Fermion-gluon vertices are also identical to those in the A^ 
fiavor Wilson fermion. Denoting the interaction vertex with n gluon as V^*-"^ we write down 
three of them here, 



V^^\p,k) 



11, st 



-igT" 



{-p + k)f, . . {-p + k)f, 
7m cos ir sm 



6 



st; 



(11.27) 






2! 
3 iT^T'^T' 



Z7^ sm — r cos 



Ofii/Ost, 



Isym. 



3! 



i-P + k)^ . . i-p + k), 
7^ cos — tr sm 



(11.28) 
S,uSf,x6st, (11.29) 



where k and p represent incoming momentum through the fermion hne. T° is a generator of 

jsym. 



color SU{Nc) and {T°'T^T'^) sym. represents summation over symmetric order of indices a, 6, c. 



Among several gluon interaction vertices from the pure gauge part we need an explicit 
form of the self interaction vertex of three gluons A" A\,, A\ with incoming momentum k, /, p 
respectively, 

G^^\k, l,p)^,, = 2igr''^!^5,, sin ^{p - l)^ cos ^ 

+5x^ sin -(A; - p)^ cos ^ + (5^^ sin -(/ - A;) a cos ^ |, (11.30) 

where /'^'"^ is the structure constant of SU{Nc) algebra. G^^^ is an odd function in gluon 
momentum. 



III. QUARK SELF-ENERGY 

We consider the loop correction to the quark propagator given by a diagram in Fig. |l|, 
where blob represents a quantum correction of some loop level. This Green function is 
written in terms of the IPE vertex function times external quark propagator with help of 



{q{p)q{-p))i 



Wn 



+ 1 



2n 9 

WqJTTI 



T.g{p,m) 



Wn 



+ 1 



WQ)m 



(III.l) 



The quark self energy Sq(p, m) can be expanded in power of the external momentum p and 
mass m keeping the logarithmic dependence on them in coefficients, 



Sg(p, m) = So + p^Si^ + mEi + p^S2 + m^S2 + mp^E2^ + 0(0^ 



(111.2) 



where the coefficients Sq, Si^, Si, S2, S2, £2^ are functions of log \p\, log \m\ and p^/m, 

(III.3) 



Si = Si(log|p|,log|m|,p/m). 



Here p^/m is introduced for generality. In ( [111.2| ) the first term represents the additive 
mass correction, the second and third term contribute to the quark wave function and 
multiplicative mass renormalization factors. The 0{a) errors arise from the next three terms. 
The omitted terms give higher order errors of (9(a^). The fact that Sq, Si^, Si, S2, S2 and 
S2^ are even functions of [p, m) is important for our discussion. Our viewpoint is that if 
Sg(p, m) is an odd function of {p, m) the additive mass term and 0{a) terms do not appear 
from the first in the expansion ( |111.2D and only the terms proportional to p or m survive to 
contribute to the corrections of wave function and mass multiplicatively. In this case, the 
leading discretization errors are (9(a^). 



In the following, we represent the self-energy in the form of loop integral as 



where /; ...^ is the shorthand of n loop integral 

r- d% d\ 



i,...U J-n (27r)4 {27ry 



(III.5) 



which is invariant against the sign flip of ^i, ■ ■ ■ ,/„. We concentrate to show the oddness 
of cr(/i, ■ ■ ■ , /„;p, m) for the variables /i, ■ ■ ■ , /„,P, "o^- If it is vahd, Sq(p, m) becomes odd 
function of {p,m): 

T,q{-p, -m) = / cr(/i, ■■■,/„; -P, -m) 

a{-li,---,-l„-,-p,-Tn) 



= -S,(p,m). (III.6) 

We begin by considering the one loop contributions, which are given by "tadpole" and 
"half-circle" diagrams of Fig. ^ (a) and (b) respectively. As was discussed in the above 
a{l]p,m) is given by multiplying the self energy of domain-wall fermion by the external 
quark line factors [C,eLs —C,oRs], [Rs^e — Lg^o] of the propagator ( [1I.23| ) ( |1I.24D . Each diagram 
gives the following integrand, 

a*^^P°i^(/;p,m) = [ee^-^,i?,](p,m)02)(_p^p)^^^^^[^^^^_^^^j(p^^)^^^(;)^ (III.7) 

{l;p,m) = l^eLs - ^oRs]{p,m)V^^\-p,l)f,^stSFil,m)tuV'^^\-l,p)i,^uvG^u{p - I) 

X [R,^,-m{p,m). (111.8) 



a 



„half— circle 



We take summation over flavor index step by step from the left paying attention to the 
external quark line factor [^eLs — ^oRs] (p, m) . A characteristic property of the factor is that 
an even and odd function is separated by the damping factor L{p)s and R{p)s- We can see 
that this separation property is preserved under each steps in the diagram. 
By using the commutation relation 

L{p)sl^ = l,R{p)s (III.9) 



and the form of internal fermion propagator in ( |II.8| ), the multiplication of the damping 
factor with the fermion propagator becomes, 

L{p)sSf{1 m)st = feL{p)t + KL{l)t + foR{p)t + KR{1\, (III.IO) 

R{p)sSf{1, m),t = feR{p)t + h'M)t + foL{p)t + KL{l)u (ULll) 

where 



fo = ilti sin lf,X, 

/e = (Wil) - e°(P))X, 



K = ^7m sin l^iY2 + 
_ im-ff, sin I ^ 



m{W{l) - e°(p)) ^e"(p)-"(') 



F(/)(e°«_e-°(p)) 



+ 



' F(/)(e°(0-e-°(p)) 



+ (vr(/) - e"(p))ri + 






m 



h'^ = i'j^ sin l^Yi + 



m{W{l) - e-"(p)) m 



h' 



im-f^ sin /^ 



F(Z)(e"«-e-"(p)) 



F(/)(e"W-e-"(p)) F(/)' 






."(0 



-a(l)^ 



Y, = A{1) 
Y, = Ail) 



;i-m2)(l-M/e-"W) 1 

F(/)(l - e-"(p)-"(^)) 1 _ e-«{p)+a{0 
:i-m2)(l-iye"«)e-2"« i 



F(0(l-e— W-°W) 



1 



-a{p)+a{l) 



III.12) 
III.13) 

III.14) 

III.15) 
III.16) 
III.17) 

III.18) 
III.19) 
III.20) 

III.21) 



Among several coefficients in ( [[II.10| ) and (pi.ll| ), /e, hg and h'^ are even functions and 
fo, ho and h'o are odd. Multiplication of the external line factor with the internal quark line 
becomes 



[^,Ls - ^oRs]ip)SFil,m)st = iUe - Uo)L{p)t + {iehe ' ^oK)L{l)t 

+ iUo - Ue)R{v)s + {^eho - ioK)R{l)t. 



^111.22) 



It should be noticed that the even-odd structure in which L^'s have even coefficients and 
i?s's have odd ones is preserved. 

Multiplication with the interaction vertex V^'^^ is given by 



[ieL,-ioRs]{p)V^'\p,k)^^,t = 
[ieLs-ioRs]{p)V^^\p,k)^,^,t 



-ig{Uf,oLt + u^eRt)ip), 
g^{Uf,eLt + u^oRt){p)S^u, 



l^eLs - ^oRs] {p)V^^\p, k)f,^x,st = ig^{u^,oLt + U^eRt){.P)SiMuSiMX, 



where u^e and m^o are even and odd function, 



u^e{p, k, m) = 7^^e(p, m) cos ■ 
u^o{p, k, m) = -ff,^o{p, m) cos ■ 



— + ^rt,o{p, m) sm ■ 



-P + k)^ 



-P + k)^ {-p + k)f, 
+ irt,e{p, m) sm -. 



(III.23) 
(111.24) 
(III.25) 



(III.26) 
(III.27) 



and we omitted the color factors. In general, the interaction vertex with even number of 
gluons \/(^") preserves the structure of damping factor and that with odd gluons V^(2n+i) 
flips it so as to assigning odd coefficients to L^'s and even to i^^'s. 

After several multiplication in the diagram the external line factor from the left meets 
the other factor from the right and flavor index is summed over using the formula. 



9 



L{pi)sL{p2)s = R{j)l)sR{j>2)s = ^ _ g,^(p,)„^(p,) = C{pi,P2), 

L{pi)sR{p2)s = 0. (III.28) 

We apply this scenario cr*^*^P°^'^(/; p, -m). The left factor [^e-^s — ^oRs\{PiT^) is multiplied 
to V^'^^ but this does not change the even-odd structure according to ( [111.24]) and then meets 
the right factor and becomes an odd function, 

[CLs - CRs] (P, m) [Rti, - UQ (p, m) = -C{p, p) {i% + ^Q . (III.29) 

In o-haif-circiej^^. ^^ ^^ -(.j^g gjp Qf ^];^g eveu-odd structure occurs twice at each interaction vertex 
V^^'^ with ( |111.23| ). Multiplication with the fermion propagator does not change the structure 
as was shown in (|TL22|). Finally the left factor meets the right factor keeping the original 
structure and gives an odd function, 

[CLs - CRs] (p, m) [Rtie - UQ {p, m) = -Cip, p) [Cio + Cie) . (111.30) 

As a result, it is shown that o"(/; p, m) is an odd function of (/, p, m) at one loop level together 
with the fact that gluon propagator is an even function. Consequently the additive mass 
correction and 0{a) errors vanish. 

It is straightforward to extend the above discussion to two loop level. As examples we 
consider the diagrams of Fig. ^ (a), (b), whose integrands are 

a^"'\li,l2\p,m) = [^eLs-^oRs]{p,m) 

xV^'\-p,h)^SFih)sM'\-h,l2)pSFil2,m)tuV^'\~k,h).SFih,m)u. 

xV^'\-lup),G^,{p-h)Gp^{h-l2)[Rv^e-L,mp,m), (III.31) 

a^^\lij2;P,m) = l^^Ls -^oRs]ip,m) 

xV^'\-p,l2)^SFil2,m)stV^'\-l2ji)uSFih,m)tuV^'\-h,p)p 

X G^'^\-p + h^h - h^P - ll)af3'yG ^a{p - l2)Gyj3{li -l2)Gp^{p-h) 

X [R,^,-L,^{p,m). (III.32) 

There are four interaction vertices V^-^^ in a^"-* (/i, I2; p, m). Since the fermion propagator does 
not affect the even-odd argument, the even-odd structure of the left factor [^e-^s— Co-Rs](p, rn) 
is flipped four times and finally meets with the right factor in the same form as in the tree 
level 

[i',L,-CR,][R,i,-L,Q (III.33) 

and gives an odd contribution. Combined with the gluon propagator that is an even function, 
we can see a^^\li,l2',p,'m) is odd. On the other hand, there are three interaction vertices 
V^^^ in a^^\li,l2]p,m). This fiips the even-odd structure of the left factor three times and 
the external quark line factor gives even contribution, 

[CLs~CRs][Rv^e-L,U (ni.34) 

However, diagram a^^\li^l2]p.,m) also contains the three gluon self interaction vertex G^^^ 
which is an odd function of (/i, hiP)- Multiplying this term, the total loop integrand is an odd 

10 



function. Confirmation of oddness for other remaining two loop diagrams is straightforward 
with the help of ( |1II.22|) - (PI.25| ). We leave it to reader and just pick up three of them in 
Fig. I (c), (d) and (e) as examples. The above procedure is applicable to any quark self- 
energy diagram at any loop level and shows the oddness of the loop correction. The 0{a) 
errors vanish for any diagrams. 

In the end we have two comments on two loop diagrams Fig. ^ (f), (g) with fermion loop 
in the gluon polarization. The loop correction to the gluon polarization does not spoil the 
evenness of the gluon propagator because it is protected by gauge symmetry. Although it 
is almost trivial for the DWQCD from the same reason, we checked it explicitly at one loop 
level. This is accomplished by evaluating 



V(Z,p) = -tr \V^'\-1, 1 + p)^Sf{1 + p, m),tV^'\-l - p, 1),Sf{1, m), 

V^^\-lJ)^,SF{l,m)J. (III.35) 



-tr 



Since the evenness is seen easily by substituting (|II.8| ), (|II.27|) and ( pi.28|) , we omit the 



explicit calculation. The oddness of the three gluon self interaction vertex is also protected 
by gauge symmetry. The second comment is concerning the Pauli-Villars field needed to 
settle the infra-red divergence in fermion loop with infinite number of domain-wall fermion. 
The Pauli-Villars field is introduced as the A^ fiavor Wilson Dirac boson whose Dirac operator 
is the same as that of the domain- wall fermion except for the physical quark mass is changed 
to the cut-off order and opposite signature; m ^ —1 [^ and it keeps the gauge invariance. 
This does not change the even-odd structure of the Feynman rules in pure gauge part. 

IV. QUARK BILINEAR OPERATOR 

We consider quark bilinear operators in the following form 

Or{x) = q{x)Tq{x), F = 1, 75, 7^, 7^75, cr^^. (IV.l) 

We calculate the quantum correction to the Green's function {Or{x)q{y)q{z)) with external 
momenta p,p', As in the previous section we see that the external line is essentially written 
in terms of the quark propagator times damping factors in ( |11.23| ) and ( [11. 241 ) . Making use 



of this fact the full Green's function for small external momentum becomes 

^0) ^wZ2 



'l-wD Z,„Zo „ , , ,„ il-wl)Z,„Z'), 



tf + (1 - WQJZ^Z^^m tf + {l- WQJZ^ZJ^m 

where Tr is the loop correction to the operator vertex and Z2, Z^, Z^ are quantum correction 
to quark wave function, quark over all factor wq and mass whose explicit value is given in 
Ref. 1^ at one loop. From the result of previous section, it is proven that the 0{a) errors 
vanish automatically for Z2, Zy^, Zm- In this section we show that for Tp. 

As in the previous section we expand the effective vertex Tr{p,p',m) in terms of ex- 
ternal quark momentum p, p' and mass m keeping the logarithmic dependence on them in 
coefficients, and then estimate the lattice artifact: 

Tr(p',p,m) =To+p^Ti^+p;T{^ + mfi + 0(a2). (IV.3) 
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The coefficients To,Ti and T/ are functions of log |p|,log |p'|, log |m|,p/m,p'/m. Tq con- 
tributes to the renormalization factor of the bilinear operator and following three terms 
are 0{a) errors. If TY{p',p,m) is an even function of {p',p,m), the 0{a) errors vanish 
automatically. In the following we adopt the loop integral form as 

W.p.m)=/ -,ik,-;L;p',P,m)- (IV.4) 

Jll---ln 

Tr is even provided the integrand ■j{li, ■ ■ ■ ,ln', p', P, fn) is even function of (/i, ■ ■ ■ , /„, p, p', m). 
At one loop level, ■j{l;p',p,m) is given by a diagram of Fig. |Kc) as 

7(/;p',p,m) = [^,Ls - ^oRsW,m)V^'\-p',p' + l)[R,r]o + L,7],]{p' + l,m)T 

X [r]oLt + r],Rt]ip + l,m)V^'\-p - l,p)[Rt^e - LtU{p,m)G^,{l), (IV.5) 



where the quark propagators ([II.14|) and ([II.15|) are used as internal lines, since "physical" 



fields are living in the operator vertex. Each internal line propagators contain damping 
factors La, Rg and consideration of even-odd structure is carried for left and right hand side 
of fermion line independently. There is one fermion interaction vertex V^^^ in each fermion 
lines of ( |1V.5| ) and flip of even-odd structure occurs once. This gives odd functions, 

[CLs-CRs][RsVo + L,r,,] (IV.6) 

for the left side line and 

[ii'M + tRtWRtie - UVo] (IV. 7) 

for right side one making the total contribution even. 

For two loop diagrams, we consider Fig. ^(a) as an example. There are two interaction 
vertices V^'^'> in each fermion line, where the even-odd flip occurs twice giving odd contri- 
bution. Therefore, the total contribution from the multiplication of them becomes even 
as expected. The above argument is valid for all other remaining two loop diagrams. For 
example estimation of diagrams in Fig. ^(b),(c) is an easy task and we leave it to reader. 

V. CONCLUSION 

In this article, we made an estimation of the lattice artifacts in loop correction perturba- 
tively for DWQCD in A^ ^ oo limit. We studied quantum corrections to the Green function 
constructed from the physical quark field only. We revealed the even-odd structure of the 
diagrams in terms of the external quark momentum and mass. 

With the help of the damping factors in the external quark propagator, we have shown 
to two loop level that the quark self-energy is definitely odd and the effective vertex of the 
quark bilinear operator is even. This means that the quark self-energy has neither additive 
mass correction nor 0{a) errors, and 0{a) errors in the effective vertex of the quark bilinear 
operator vanish automatically. Since our discussion can be extended to any loop diagram, 
we can show that 0{a) errors do not appear for any diagram at any loop level in the same 
way. This is because the disappearance of them is deeply concerned with the hidden axial 



symmetry |18| of the DWQCD which was explicitly shown to exist in Ref. [21 1 for effective 



Lagrangian of the physical quark fields. 
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FIGURES 




FIG. 1. A quantum correction to the quark self-energy at some loop level . The blob represents 
the IFF part. 
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P v(2) P 



V<i> Sp(l) yd) 



(a) (b) (c) 

FIG. 2. The one loop contributions to the quark self-energy; (a), (b) and the quark bilinear 
operator; (c). 






(f) (g) 

FIG. 3. Typical contributions to the quark self-energy at two loop level. 
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(a) (b) (c) 

FIG. 4. Typical contributions to the quark bilinear operator at two loop level. 
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